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Many-body systems relaxing to equilibrium can exhibit complex dynamics even if their steady 
state is trivial. At low temperatures or high densities their evolution is often dominated by steric 
hindrances affecting particle motion [1,2,3]- Local rearrangements are highly constrained, giving 
rise to collective - and often slow - relaxation.This dynamics can be difficult to analyse from first 
principles, but the essential physical ingredients are captured by idealized lattice models with so- 
called kinetic constraints [4]. Here we experimentally realize a many-body system exhibiting manifest 
kinetic constraints and measure its dynamical properties. In the cold Rydberg gas used in our 
experiments, the nature of the constraints can be tailored through the detuning of the excitation 
lasers from resonance [5,6,7,8], which controls whether the system undergoes correlated or anti¬ 
correlated dynamics. Our results confirm recent theoretical predictions [5,6], and highlight the 
analogy between the dynamics of interacting Rydberg gases and that of soft-matter systems. 

PACS numbers: 34.20.Cf, 32.80.Ee 

Complex collective relaxation in many-body systems is often accompanied by a dramatic slowdown of diffusion 
processes and the emergence of non-ergodic and glassy phases [1,2,3]. These features can be seen to be the consequence 
of effective kinetic constraints in the dynamics [9]. A kinetic constraint is a condition on the rate for a local transition 
dependent on the local environment: the transition and its reverse - irrespective of whether they are energetically 
favourable or unfavourable - can only occur if the constraint is satisfied. Kinetic constraints can severely restrict 
relaxation in situations where local particle arrangements make satisfying them unlikely, which is typical of fluid sys¬ 
tems with excluded volume interactions such as dense colloids or supercooled liquids. When a constraint is satisfied, 
however, the transition is allowed and a local rearrangement is ’’facilitated” [4,9]. Kinetic constraints naturally give 
rise to collective and spatially heterogeneous relaxation, and are used to describe situations where the correlation 
properties of the dynamics go beyond those of the static stationary state, a salient feature of glassy systems [3]. While 
kinetic constrains and dynamic facilitation are argued to play a central role in the behaviour of glass formers [9], 
it is difficult unambiguously to establish a connection between microscopic processes and emerging kinetic constraints. 

In our experimental study we report the emergence of manifest kinetic constraints in the many-body dynamics 
of laser-driven gases of Rydberg atoms in the presence of noise. We explore two kinds of constraints which lead 
to the suppression or the facilitation of atomic excitations in the vicinity of atoms excited to Rydberg states. Our 
experimental technique allows us to measure both the mean number of excitations and the fluctuations around the 
mean, and we find clear signatures of the constraints in the pronounced (anti-) correlations evident in both quantities. 
The experimental observations are in excellent agreement with numerical simulations of a kinetically constrained 
system of Ising spins, where the constraint is explicitly derived from the actual interatomic interactions. 

The kinetic constraints and resulting correlated excitation processes explored in this paper are shown schematically 
in Fig. 1. We consider a many-body system consisting of atoms that can be represented as pseudo Ising spins having 
a (non-interacting) ground state | g) (spin down) and an excited state |r) (spin up) coupled by a Rabi frequency fi. 
Atoms at positions r, and r 3 in the Rydberg state |r) interact through the van der Waals interaction V tl = j r j 6 , 
where Cq is the van der Waals interaction coefficient [10,11] (which, in the present experiment, is positive). Here 
we explore the regime of incoherent evolution, 7 >> fl, where 7 is the decay rate of the atomic coherences. In 
this strong dissipation limit the evolution of the many-body system is restricted to the subspace of classical spin 
configurations by virtue of the quantum Zeno effect [12]. The dynamics can be written as a classical rate equation 


with (de-)excitation rates R(A) = 


1+ 



t [ 5 , 7 ], where rij = 1 if the atom at position r 3 is in the 
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Rydberg state and 0 otherwise. Generally, therefore, the excitation rate I\(A) for an atom depends on the state of 
all the atoms in its vicinity and on the detuning A of the excitation laser from resonance. For A = 0 (in practice, 
A < 7 ), the interactions between an atom and its neighbours lead to a blockade constraint resulting in anti-correlated 
dynamics: the more excited atoms there are in the vicinity, the smaller T,, leading to a spatially inhomogeneous local 
excitation rate and an overall slowing down of the dynamics as the number of excitations in the systems grows. The 
inter-particle distance below which this blockade constraint becomes important is given by the (incoherent) blockade 
radius r b = (Cg/h"/)® . By contrast, for A > 0 the single-atom excitation rate is small, but now a facilitation 
constraint appears: an excited atom somewhere in the system shifts atoms contained within a shell of radius rf ac and 
width 5rf ac into resonance, where rf ac = (Ce/HA)i , which corresponds to the van der Waals interaction compen¬ 
sating the laser detuning [8,13,14,15,16,17], and Srf ac = 7 ^ 7 . This facilitation leads to strongly correlated dynamics. 

In our experiments the above model is realized with 8 'Rb atoms in a magneto-optical trap (MOT) that are excited 
to repulsively interacting 70S Rydberg states [18] and detected by a channeltron after field ionization (see Methods). 
In order to verify that the dynamics of our system is in the regime of validity of the above rate equation, we first 
measure the resonant excitation dynamics, represented by the mean number of excitations A as a function of time t 
(here and in the rest of the paper we report the derived actual quantities obtained by dividing the observed quantities 
by the detection efficiency), for different values of the Rabi frequency Cl. In the incoherent regime (7 >> Cl) a 
scaling of Tj, and hence of the overall dynamics, with Cl 2 / 7 is expected (as opposed to the coherent regime, where 
the theoretical linear scaling with Cl was demonstrated in [19]). Fig. 2 shows the results of experiments for three 
different values of Cl. When plotting the experimental data versus t ■ Cl 2 / 7 all three curves collapse onto a single 
curve, demonstrating the expected scaling of the incoherent dynamics. 

We now proceed to investigate the constrained non-equilibrium dynamics of our system. To realize the blockade 
constraint we set A = 0 and measure A as a function of time. In order not to saturate our detection system for 
long excitation times (which, in practice, means keeping the observed number of excitations < 40) we use the ID 
configuration (see Methods). The degree of correlation of our system is controlled by the number of atoms per 

blockade length (where a = 3 is the mean distance between N g ground state atoms in the excitation 

volume V exc , and 77 = 11.1 //m for our parameters), which we vary by changing the effective density of the MOT (see 
Methods). In this way, we can prepare samples with atom numbers corresponding to values of — between around 
1.3 (i.e., close to the non-interacting case ^ ^ 1 ) and ^ = 4.2. The results of these experiments are shown in Fig. 
3 a), together with a numerical simulation (see Methods) that exhibits excellent qualitative and good quantitative 
agreement (to within overall factors between 0.5 and 2, which are indicated in the captions to Figs. 3 and 4). The 
cross-over between the initial excitation regime in which the interactions play no role, reflected by a linear increase in 

A, and the blockade regime is clearly visible, indicating the point at which the average distance d = 3 between 

excitations becomes smaller than 77 [20]. When the results are normalized by the number of ground state atoms N g 
in the excitation volume, which yields the fraction of excited atoms, in the blockade regime the excited fraction levels 
off much more sharply for large values of ^. 

The above observations can be represented more quantitatively by examining the average growth rate of excitations 
per atom, ( dN/dt)/N g , as a function of d. The growth rate is extracted by numerically differentiating the A vs 
t data from Fig. 3 a) after smoothing them in order to avoid artefacts due to noisy data. We expect this rate to 
decrease sharply with decreasing values of d below the blockade radius, which defines the effective range of influence 
of a Rydberg atom. This is clearly borne out by experiment, as can be seen in Fig. 3 c), in which all the data sets 
follow a single curve which is flat for d > 77 and decreases by four orders of magnitude between 11 /jm and 6 /im. 
We compare this universal curve to the theoretically predicted average single-atom growth rate derived from the rate 
function Tj(A = 0) in a mean-field approach [5] and find good agreement. 

The strong correlations caused by the kinetic constraints are expected to also be clearly visible in the fluctuations 
of the systems, which we investigate through the Mandel Q-parameter [21]. Similarly to the dipole blockade in the 
coherent regime [22,23], the effect of the blockade constraint is expected to reflect itself in a negative Q-parameter 
that decreases with the number of excitations since for large numbers the system has fewer choices for distributing the 
excitations, which in turn leads to reduced fluctuations. Moreover, following the above reasoning, the Q-parameter 
should depend exclusively on the number of excitations and not on the number of ground state atoms, so plotting Q 
as a function of A for all three data sets should yield a single curve. Again, this is experimentally confirmed in Fig. 3 d). 
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In contrast to the blockade constraint, which causes strong anti-correlations in the dynamics, the facilitation con¬ 
straint should lead to a strongly correlated evolution. In order to explore this regime, we now set A/27T = +19 MHz, 
for which Tf ac = 6.4/zm and 5rf ac = 39 nm. Since we expect the predicted facilitation dynamics to be the more 
pronounced the larger the overall facilitation volume, which grows with an increasing number of excitations, we 
choose the 3D configuration for this experiment. The results are shown in Fig. 4 (a). In this figure, three stages 
can clearly be distinguished: the initial nucleation stage (for t < 10 /zs) in which N grows slowly due to off-resonant 
single particle excitations that are suppressed by a factor « 1.4x 10~ 3 compared to the resonant regime, 

and in most cases no seed excitation has yet been created to nucleate the facilitation process; the facilitation stage 
(10 /zs < t < 50 /zs) in which the number of excitations grows fast due to successive facilitation starting from an 
initial seed; and a saturation stage (t > 50/zs), in which the dynamics decelerates again due to the finite size of the 
atomic cloud (this regime is visible in the experimental data but already affected by spontaneous decay, which is 
not included in the simulations; the lifetime of the 70S state is around 150/zs). These regimes are also reflected in a 
plot of the average total growth rate as a function of time (Fig. 4(b)), in which the acceleration of the dynamics at 
around 10 /zs can be clearly seen. 

To demonstrate that the sign of the detuning is crucially important, we repeated the above experiment 
for A/27T = —19MHz (Fig. 4(a); for comparison with the blockaded constraint, data for A = 0 are also shown). 
Here, the facilitation condition cannot be fulfilled, and only off-resonant excitation of single Rydberg atoms is possible. 

Apart from the acceleration of the dynamics in the facilitation stage, the strong correlations induced by the 
facilitation constraint are clearly discernible in the Mandel Q-parameter shown in Fig. 4 (c). For A/27T = +19 MHz, 
Q grows up to 30 /zs, becoming large and positive, and then decreases to around 0. This can be understood as 
follows: in the facilitation regime an excited atom enables the probabilistic excitation of further atoms, so that small 
fluctuations in the timing of the excitation are amplified for short times, leading to a large value of Q. For long times, 
the dynamics becomes similar to the blockaded case as no further facilitation is possible, and hence Q decreases. This 
non-trivial behaviour of the correlations is a genuine many-body effect and, therefore, absent for negative detuning. 
The fact that Q is not 0 in that case, as it should be for a purely Poissonian process, is probably due to technical 
noise and fluctuations in the atom number, laser intensity and detuning. 

We have shown that manifest kinetic constraints govern the dynamics of Rydberg gases. The constraints are a 
consequence of quantum mechanics, but in the strong dissipation regime studied here the dynamics is effectively 
classical. Future experiments could systematically probe the crossover from the incoherent to the coherent regime, 
giving insight into the role of quantum effects in constrained dynamic relaxation. Furthermore, techniques for 
spatially resolving Rydberg excitations [24,25,26] should reveal rich spatial correlations in the dynamics [27], in 
analogy with dynamic heterogeneity in glasses [3]. More generally, our results give a first indication of the broad 
potential for applying ideas from soft-matter physics to the study of interacting atomic systems. 
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METHODS 
Experimental protocol 

We create cold samples of 8 ^Rb atoms in a magneto-optical trap. The resulting roughly spherical atomic clouds 
have a Gaussian density profile with widths between 45 and 160/zm and contain between 10 4 and 1.2 x 10 6 atoms. 
We control the (effective) number of atoms and hence the density of the clouds (keeping their size and shape 
constant) by depumping a fraction of the atoms in the F=2 hyperfine ground state to the F =1 state using a resonant 
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laser pulse of duration up to 2ps (for details see [28]). This state lies 6.8 GHz below the F =2 state (the ground 
state for the Rydberg excitation process) and hence does not couple to the Rydberg excitation lasers. After the 
depumping pulse, we infer the effective atom number by exciting a few Rydberg atoms (less than 10, in order to 
avoid interaction effects) on resonance and measuring the total growth rate, which in the non-interacting, incoherent 
regime is proportional to the atom number. 

Excitation to the 70S Rydberg state (for which Cq = hx 1337 GHz pm 6 ) is achieved by a two-photon process with 
laser beams near 420nrn (waist 10 or 40pm) and 1013 nrn (waist 110pm), and the 420nrn laser detuned by 660MHz 
from the 6 P 3/2 excited state. Depending on the sizes and propagation directions of those beams, we can realize two 
geometries: and effective ID geometry, in which the 420 nm laser is focused to a waist of around 10 pm (comparable 
to the blockade and facilitation radii) and the 1013 nm intersects the 420 nm beam and the atomic cloud at a 45 
degree angle; and a 3D geometry, in which the 420 nm laser has a waist of 40 pm and co-propagates with the 1013 nm 
laser. The combined decoherence rate due to the laser linewidths and residual Doppler broadening is estimated to be 
7 = 27r x 0.7 MHz. 


After an excitation pulse of up to 100 ps, the Rydberg atoms are field ionized and detected by a channeltron with 
an overall detection efficiency of around 40%. The experiments are repeated 100 times for each set of parameters in 
order to obtain the mean number and variance. 


Numerical simulation 

The dynamics of a laser driven Rydberg gas in the limit of strong dissipation is given by the equation [12] 
v = r*(l r »)(5i| + \gi){fi\ — \ri){rt\ — \gi)(gi\)v where v encodes the probability of the system to be in one of its 

classical configurations (products of the local spin states |r;) and | gij). In order to account for the inhomogeneity 
of the atomic cloud we distributed the atomic positions according to a Gaussian profile in close approximation to 
the experimental situation. The spatially inhomogeneous laser profile was accounted for by using a local varying 
Rabi frequency for the calculation of the value of the rate function Tj for each individual atom. In the simulation 
only atoms inside an effective excitation volume V exc were considered. This excitation volume was defined as the 
region in which the free flipping rate of an atom was at least 10% (3D data) or 5% (ID data) of the maximum rate. 
The simulations where conducted by implementing a classical kinetic Monte Carlo algorithm, and the results were 
averaged over 1000 runs (ID data) or 100 runs (3D data). 


Mean-field calculation of the excitation rate 

The theoretical curve in Fig. 3 d) was obtained from the expression for Tj(A) in a mean-field approach by averaging 
fi 2 over the excitation volume defined above and substituting |r» — r : j | by half the mean distance between excitations, 
where the volume used for calculating the mean distance was chosen as 1.92 x 10 4 pm 3 , which gave the best agreement 
with experiment. The interaction term in Tj(A) was multiplied by 6 in order to account for six nearest neighbours 
for each atom. 
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Figure 1 Kinetic constraints realised in a gas of Rydberg atoms, a) Schematic representation of the laser-induced 
coupling between ground and Rydberg states with Rabi frequency Q and detuning A, where coherence between the 
atomic states is lost at a rate 7 . Each atom can be effectively described as an Ising pseudo spin. In the incoherent 
regime 7 >> f2, the dynamics of an atom reduces to incoherent state changes (spin flips) at a rate Tj(A). b) 
Interaction potential V(| 77 — 77 |) between two excited atoms at positions 77 and 77 (black line) and excitation rates 
r^(A) for A > 0 (blue line) and A = 0 (grey line). The excitation rate for A = 0 drops to zero for interatomic 
distances 1 77 — 77 1 below the blockade radius 77 ( blockade constraint ), whereas it peaks at the facilitation radius rf ac 
for A > 0 (facilitation constraint). This leads to blockaded dynamics, shown in c), and to facilitated dynamics, 
shown in d), respectively. In c), excitation of individual atoms occur at the resonant rate until the distance between 
adjacent excitations approaches the blockade radius. The third spin from the left is excited at a strongly reduced 
rate as it is within the blockade radius. In d), a seed excitation enables further excitations at the facilitation distance 
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Figure 2 Dependence of incoherent dynamics on the Rabi frequency. The inset shows the mean number of excitations 
in a sample of 7 x 10 5 atoms in a MOT of width 60 /rm in the 3D configuration as a function of time for three different 
Rabi frequencies: D/27T = 81 (open diamonds), 43 (green squares) and 20 kHz (red circles). When multiplying 
the excitation times by SI 2 /7 (where 'y/2n = 0.7MHz), the three curves collapse onto each other (main figure), 
demonstrating the expected Sd 2 scaling in the incoherent excitation regime. Representative error bars are one s.d.m. 
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Figure 3 Evidence for the blockade constraint in a gas of Rydberg atoms, a) Mean number of excitations and b) 
excitation fraction as a function of time for different atom numbers N g inside the interaction volume: 5600 (open 
diamonds), 715 (green squares) and 180 (red circles). These correspond to values of the parameter ^ of 4.2, 2.1 
and 1.3, respectively. The numerical simulations (solid lines) have been scaled vertically by a factor of 1.8. c) 
Average single-particle excitation rate as a function of the mean distance d between excited atoms. The dashed 
line is obtained from the expression for Tj(A) in a mean-field approach (see Methods), d) Mandel Q-parameter as 
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a function of mean excitation number. The vertical axis on the left corresponds to experimental data, that on the 
right to the numerical simulation (solid lines). The deviations of the experimental data from the simulation are most 
likely due to experimental noise (positive Q for small values of N) and saturation effects of the detection system 
(values of Q < — 1 for large N). Representative error bars are one s.d.m. 



Figure 4 Evidence for the facilitation constraint in a gas of Rydberg atoms, a) Mean number of excitations as a 
function of time for A/27T = +19MHz (blue circles), A/27T = —19MHz (red circles) and A = 0 (grey circles). The 
numerical simulations (solid lines) have been scaled vertically by a factor 0.56. Representative error bars are one 
s.d.m. b) The growth rate of excitations as a function of time. The left-hand vertical axis refers to the data for 
A/27T = +19 MHz and A/27T = —19 MHz, the left-hand axis refers to A = 0 . The numerical simulations (solid 
lines) have been scaled vertically by a factor 0.56. c) The Mandel Q-parameter as a function of time. The numerical 
simulations (solid lines) have been scaled by a factor 2. 
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